Probing Quark-Gluon Interactions with Transverse Polarized Scattering 
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We have extracted QCD matrix elements from our data on doubly polarized inelastic scattering 
of electrons on nuclei. We find the higher twist matrix element d2, which arises strictly from quark- 
gluon interactions, to be unambiguously non zero. The data also reveal an isospin dependence of 
higher twist effects if we assume that the Burkhardt-Cottingham Sum rule is valid. The fundamental 
Bjorken sum rule obtained from the ao matrix element is satisfied at our low momentum transfer. 



PACS numbers: 11.55.Hx,25.30.Bf,29.25.Pj,29.27.Hj 

Polarized deep inelastic scattering (DIS) of leptons on 
nucleons has helped to establish the modern understand- 
ing of nucleon structure. Measurements of the longitu- 
dinal nucleon spin structure function (SSF) gi in inclu- 
sive doubly polarized DIS, where only the final lepton 
is detected, have revealed that the nucleon spin does 
not result from the naive sum of the spins of its three 
valence quarks. Rather, contributions from sea quarks, 
gluon spin and orbital angular momentum seem to be 
required [H-Q. DIS is well understood at high energies, 
when the virtual photons exchanged between the incident 
beam and the target illuminate the nucleon constituents 
as the asymptotically free quarks of QCD. This is the 
regime of small separations between quarks and giuons, 
which is best tested with short wavelength electromag- 
netic probes. But the processes that polarized DIS ex- 
plores are dominated by leading order terms in the expan- 
sion of gi that are related to the non-interacting partonic 



description of QCD. For a full understanding, the inter- 
actions among nucleon constituents must also be probed 
by looking into contributions beyond leading order, which 
grow in importance at lower energies. Of particular inter- 
est are the so called 'higher twist' terms, which measure 
correlations among partons (quarks and giuons), and of- 
fer unique perspectives on the phenomenon of confine- 
ment inside nucleons and other hadrons. 

Access to these subleading interactions is possible with 
transverse polarized inelastic scattering. In this configu- 
ration, the target spin is held orthogonal to the incident 
lepton spin, and the 52 spin structure function dominates 
compared to 171 Q ■ QCD matrix elements that arise from 
multi-parton correlations can be extracted directly from 
moments of g2. However, this approach has been lim- 
ited due to the technical challenges involved in aligning 
proton spins at 90° relative to the beam. As we detail 
below, we have carried out a complete measurement of 
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the longitudinal and transverse inclusive SSF's gi and (72 
of the proton and deuteron, with the goal of extracting 
matrix elements related to quark-gluon, qg, interactions. 

The Operator Product Expansion (OPE) fs'] relates the 
moments of the SSF's to reduced quark and giuon ma- 
trix elements representing the multi-parton correlations 
that lead to confinement. This perturbative approach 
to QCD is formulated in terms of the Cornwall-Norton 
(CN) moments Q: 



(1) 



where x = /2Mv is the Bjorken scaling variable, 
is the four-momentum transfer squared of the scattering 
process, M the nucleon mass, and v the energy transfer. 
Conventionally, the index n = 1 is not displayed. 

In the OPE, the CN moments are expanded in a power 
series in l/Q'-''"^' of increasing 'twist' t, which is defined 
as the mass dimension minus the spin n of the corre- 
sponding QCD operator The leading twist-2 terms 
map to the successful predictions of the parton model, 
such as the Q^-independence of the structure functions 
up to logarithmic corrections. Higher twist (HT) correc- 
tions arise from the non-perturbative multiparton inter- 
actions, whose contributions at low energy increase as 
l/Q'^, reflecting confinement. Specifically, twist-3 repre- 
sents qgq correlations, so a non-zero result at a given 
for the term representing twist-3 signals a departure from 
the non-interacting partonic regime. 

Studies of higher twist effects [7l-[l]| have typically fo- 
cused on the CN moments I{Q^) = Jq dx x^ {2gi + 8512) 
and ri((5^), in order to extract the twist-3 and twist-4 
matrix elements. However, as recently stressed in (6l.[l^. 
this approach is only appropriate when terms of purely 
kinematic origin, known as 'target mass corrections' 
(TMC) [13l4la | for their connection to the finite nu- 
cleon mass, can be neglected. These terms, of order 
0{M^ /Q'^), are formally related to twist-2 operators and 
must be cleanly separated from the desired dynamical 
higher twists. I{Q^) is often labeled 'o?2' in the litera- 
ture, although, as discussed below, it is not equivalent to 
the HT matrix element indicated by the same symbol. 

As emphasized in [l^], dynamical HT can be ex- 
tracted to order 0{M^/Q^) by using Nachtmann mo- 
ments [Tot , which depend on the scaling variable 
£, = 2x/{l + y/1 + {2xMy/Q'^). This variable extends 
scaling as x approaches 1, for data measured at low Q^. 
The Nachtmann moments are defined as 
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5.1 


RC 


3.2 


4.0 


9.5 


4.5 


Fi 


3.0 


4.0 


3.0 


4.0 


R 


0.9 


1.8 


1.2 


2.8 


PbPt 


1.6 


1.6 


5.2 


3.1 


dependence 


0.33 


0.1 


4.3 


0.3 


Total 


6.8 


8.0 


13.0 


8.9 



TABLE I. Measured integral systematic uncertainties (in %) 
arising from the target dilution factor /, radiative corrections 
RC, the _Fi and R data fits, beam and target polarizations Ph 
and Pt, and the evolution of our data fit to constant . 



Here, the a„ ((i„) represent twist-2 (-3) matrix elements, 
while the E2 are the corresponding Wilson coefficients, 
which contain logarithmic QCD corrections. For conve- 
nience, these corrections are absorbed into the definition 
of the effective matrix elements a„ and dn- The index 
runs over n = 1,3,... for Mi, and n = 3, 5, . . . for M2. 
When Q^^M^, the Nachtmann moments simplify: 

Consequently, deviations from unity of the ratio of Nacht- 
mann to CN moments give a quantitative measure of the 
TMCs at finite Q^: 
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The leading twist-2 contribution to 32 can be deter- 
mined from gi via the Wandzura-Wilczek relation [27} : 
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So any twist-3 effects arise entirely from g2 



„ „ww 

92-92 ■ 



Experiment EO 1-006 was conducted in Hall C of the 
Thomas Jefferson National Accelerator Facility by the 
Resonance Spin Structure (RSS) collaboration. We mea- 
sured the parallel and perpendicular double spin asym- 
metries and A± in the scattering of 100 nA po- 
larized 5.755 GeV electrons on polarized protons and 
deuterons. Scattered electrons were detected at an an- 
gle of 13.15° using the Hall C High Momentum Spec- 
trometer. The kinematic coverage in invariant mass was 
1.090 <W < 1.910 GeV, corresponding to Xq = 0.316 < 
X < Xmax = 0.823, at an average four-momentum trans- 
fer of ((52) = 1.28 ±0.21 GeV2. Systematic uncertainties 
are detailed in Table HI with more details in [20| and (2l| . 

The SSFs were extracted from the measured asymme- 
tries using a fit to the ratio of longitudinal to transverse 
cross sections R [l^l and the unpolarized structure func- 
tion Fi [l^l , based primarily on data measured previously 
in Hall C. The standard formulas and procedure used to 
obtain A|| , Aj _ and the SSFs from the data are detailed 
in [H] and [20l| . Simultaneous determination of both SSFs 
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FIG. 1. Left-top: RSS proton (72 data [20|], along with the 
RSS fit [2^. The shaded curve is g""' evaluated from the 
RSS fit. Left-bottom: RSS deuteron g2 data, and similar 
curves as above. Right-top: The open circle is the mea- 
sured RSS proton V^{Q^) data. The full circle is Ar2(Q^). 
The inner (outer) error bars represent statistical (total) un- 
certainty. The shaded curve is the elastic contribution to 
ViiQ^). Right-bottom: Same as above, but for deuteron. 



allowed us to evaluate the moments without model input 
for g2 , as was necessary in some previous analyses [9!] . 

The moments reported herein are evaluated at (Q^) — 
1.28 GeV^, and have been decomposed into contributions 
from the measured resonance region xq < x < Xthr (la- 
beled with 'm'), the well known x — 1 elastic ('el') contri- 
bution, and the unmeasured ('u') portion below xq. We 
note that the small difference between our experimental 
Xmax and the nucleon inelastic threshold xthr has neg- 
ligible impact on the integrated results. We have used 
fits [2^ to our data to evaluate the moments in the mea- 
sured region. The nucleon elastic contribution was cal- 
culated using the form factor parameterizations of p3l |. 
Relative uncertainties of 5, 1, 14, and 2.5% have been 
assumed for the electromagnetic form factors ^ and 
^EM respectively. For the deuteron, x = 1 represents 
quasielastic scattering, the strength of which we estimate 
by combining the nucleon elastic contributions using the 
D-state correction discussed below. The deuteron nu- 
clear elastic contribution is negligible here. We have also 
evaluated the neutron and non-singlet (proton - neutron) 
moments by using the relation [24] F" = F^ — F*', where 
the singlet F^ = T'^/lD, and 7^ ^ 0.926 ± 0.016 [H] 
is the D-state correction to the deuteron wave func- 
tion. The uncertainty arising from this approach is es- 
timated [il] to be 0(1%). The singlet and non-singlet 
results assume negligible heavy quark contributions. 



Table |TT] provides numerical values for the moments. 
The ratio R2 differs significantly from unity, indicat- 
ing large unwanted twist- 2 kinematical contributions to 
I{Q^)- The full results for the matrix element d2 indicate 
clear twist-3 at more than 6a for the proton and 3(t for 
the neutron. These non- vanishing results unambiguously 
indicate the presence of qgq correlations. Their magni- 
tudes can be used in comparing with nucleon models. 

The values of I{Q^) and 6,2 for x < xq were estimated 
to be zero, with a systematic error described below. In 
this unmeasured region, ^ « x and the CN and Nacht- 
mann truncated moments converge, so there is little dif- 
ference between I{Q^) and d2- I{Q'^) is expected to be 
small in this unmeasured region, because of the decreas- 
ing importance of higher twists at low a;, and the strong 
suppression due to the x^ weighting of the integral. This 
assumed low x behavior is supported by our data. Fig- 
ure[T](left) displays the x— dependence of the measured g2 
structure function. It is clear that 52 is nearly constant 
and consistent with zero within errors for the proton, 
near the low end of our measured range < a; < 0.37. 
The deuteron data shows a similar behavior, different 
from zero only at the one a (statistical) level. Neutron 
data [2^ at a; ^ 0.2 and similar are also consis- 
tent with ^2—0. We take the error on g2, evaluated 
at xo, Sg2, as a conservative upper limit for the inte- 
grand of I{Q^) in the unmeasured low x region. For the 
proton, we determined this upper limit by assuming a 
constant extrapolation of the value of 5^2(2^0) down to 
a; = 0. For the deuteron, we evaluated both a linear and 
constant extrapolation, averaged both assumptions and 
added quadratically one-half their difference as 'model' 
error to the <5^2 fit error to obtain the value in Table |lll 
A divergence of 32 as a: could invalidate this as- 
sumption. Such a possible divergence for x < 0.001 was 
discussed in [2^. Normalizing the low x dependence 01^2 
given in p9j to our 5'g2 estimates, we find the additional 
contribution to be negligible. 

It is instructive to compare our twist-3 results to pre- 
vious measurements. SLAG E155 {l\ reported an eval- 
uation of I[Q^) at (Q^) = 5 GeV^. We have corrected 
those results for TMGs [13, to obtain 0.0028 ± 0.0015 
and 0.0072 ± 0.0044 for the proton and neutron, respec- 
tively. For direct comparison, we performed a pQGD evo- 
lution 30] from RSS to the SLAG kinematics. At LO, 
we find dl = 0.0021 ± 0.0006 and 4' = 0.0031 ± 0.0038, 
which are consistent with the E155 results. The elastic 
contribution [2^ at these kinematics is smaller than the 
uncertainties and is not included in the results of this 
paragraph. NLO corrections 31] to our data have been 
calculated and will be discussed in a future publication. 

Twist-3 effects also manifest in the first moment of g2- 
The Burkhardt-Gottingham (BG) sum rule ^32,] predicts 
that F2 vanishes for all . This sum rule can be derived 
from the unsubtracted dispersion relation for the virtual- 
virtual Gompton scattering amplitudes, in an analogous 
fashion to the more famous GDH and Bjorken sum rules. 
It provides a remarkably clean tool to investigate nucleon 
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Proton 




Deuteron 


Neutron 


Non-Singlet 


a;-range 




Value 


Stat 


Syst 


Value 


Stat Syst 


Value 


Total 


Value 


Total 




Ml 


0.0676 




0.0069 


0.0274 


0.0104 


-0.0381 


0.0132 


0.1057 


0.0178 


< a; < xo 


Ti 


0.0683 




0.0069 


0.0280 


0.0099 


-0.0381 


0.0127 


0.1065 


0.0176 


(unmeasured) 


da 


0.0 




0.0008 


0.0 


0.0013 


0.0 


0.0016 


0.0 


0.0021 




/ 


0.0 




0.0008 


0.0 


0.0013 


0.0 


0.0016 


0.0 


0.0021 




Ml 


0.0330 


0.0005 


0.0022 


0.0290 


0.0010 0.0023 


-0.0016 


0.0036 


0.0346 


0.0054 


Xo < X < Xthr 


Ti 


0.0351 


0.0005 


0.0024 


0.0315 


0.0011 0.0025 


-0.0010 


0.0039 


0.0361 


0.0058 


(measured) 


d2 


0.0037 


0.0004 


0.0005 


0.0048 


r\ i~\r\i~\o r\ r\r\r\ a 

0.0008 0.0004 


0.0015 


0.0012 


0.0022 


0.0016 




I 


0.0057 


0.0006 


0.0007 


0.0082 


0.0013 0.0007 


0.0031 


0.0019 


0.0026 


0.0026 




Ml 


0.0287 




0.0006 


0.0338 


0.0012 


0.0078 


0.0005 


0.0208 


0.0008 


X ^1 


Ti 


0.0351 




0.0010 


0.0373 


0.0015 


0.0051 


r\ r\r\r\'~f 

0.0007 


0.0300 


0.0012 


(elastic) 


d2 


0.0067 




0.0008 


-0.0021 


0.0010 


-0.0090 


0.0007 


0.0157 


0.0011 




I 


0.0306 




0.0025 


0.0088 


0.0029 


-0.0212 


0.0019 


0.0518 


0.0031 




Ml 


0.1293 


0.0005 


0.0073 


0.0902 


0.0010 0.0108 


-0.0318 


0.0137 


0.1611 


0.0187 


< x < 1 


Ti 


0.1385 


0.0005 


0.0074 


0.0968 


0.0011 0.0103 


-0.0340 


0.0134 


0.1725 


0.0185 


(Total) 


d2 


0.0104 


0.0004 


0.0013 


0.0027 


0.0008 0.0017 


-0.0075 


0.0021 


0.0179 


0.0029 




I 


0.0364 


0.0006 


0.0027 


0.0170 


0.0013 0.0032 


-0.0180 


0.0031 


0.0544 


0.0045 


< x < 1 


Ri 


0.933 


0.003 


0.932 


0.012 


0.936 


0.037 


0.934 


0.008 




R2 


0.286 


0.014 


0.159 


0.078 


0.415 


0.045 


0.329 


0.025 



TABLE II. RSS Moments evaluated at (Q^) = 1.28 GeV^. The ratio i?i = MllVx and i?2 = da//. 



structure, since it is free from both QCD radiative correc- 
tions and TMCs [s^]- Initial measurements by the El 55 
collaboration [l| found a 2a inconsistency with the pro- 
ton BC sum rule at large Q^, while the same group found 
the deuteron sum rule to hold. Later measurements 0,0] 
at lower momentum transfer found agreement with both 
the neutron and '^He BC sum rules. 

The leading twist satisfies the BC sum rule ex- 
actly, so any violation must arise from higher twist ef- 
fects. For convenience we split the full integral as: 

12— ^2 +^2+^2 

= r2+r" + r^' (7) 

where the overbar signifies removal of the leading twist 
contribution, and we have made use of the fact that 
pww = Q definition. The elastic contribution is well 
known, as displayed in Fig. [TJ and we can evaluate r2 
directly from our own data. is pure higher twist, and 
cannot be directly determined, due to the lack of low x 
gi data. However, the difference, 

Ar2 = F2-r2 (8) 

depends only on measured data. A significant non-zero 
result for this quantity would indicate a higher twist con- 
tribution to the integral in the region x < xq, assuming 
of course, that the BC sum rule holds. Wc find that the 
proton results are consistent with vanishing HT effects at 
low X. On the other hand, the singlet results indicate the 
need for significant HT contributions to r2 at low a;, if 
the BC sum rule is to be satisfied. This could indicate an 
isospin dependence of the BC sum rule, or alternatively, 
a modification to the sum rule due to nuclear effects. 

Nachtmann moments are also useful for obtaining lead- 
ing twist matrix elements free of TMCs. The ap non- 
singlet matrix element corresponds to the lowest moment 



F^^" of the isovector nucleon gi, which is related to the 
nucleoli axial charge qa via the Bjorken sum rule [s^ . 
This relation is a direct consequence of QCD, and exper- 
imental tests P, [2| of the sum rule have played a criti- 
cal role in confirming QCD as the correct theory of the 
strong interaction. The target mass corrections to Fi 
given by the ratio R\ are small at our kinematics, due to 
a combination of dominant contributions from the region 
below our lowest measured xq, and the nearly identical 
values of x and ^ in that region. It is interesting to note 
that R\ « 93% for all targets, when the elastic contri- 
bution is included. The M\ non-singlet result including 
the elastic part (0.161 ± 0.019) is in excellent agreement 
with theory: Ref. [Hj predicts 0.212, which reduces to 
0.155 after application of the next to next-to-leading or- 
der corrections [111 appropriate at our Q^. For the g\ 
integrals, the contribution from the a; < region was 
calculated from Regge-inspired fits to SLAC g\ and 
data within a band = 1.3 ± q.S GeV^. We find 
gi = ax^(\ - xf{\ c/Q2) with a = 0.392 ± 0.254, 
6 = 0.0676 ± 0.084, c = 0.0636 ± 0.681 for the pro- 
ton, and a = 1.778 ± 1.971, h = 0.739 ± 0.407, and 
c = -0.786 ± 0.0256 for the deuteron. 

In summary, we have measured the spin structure func- 
tion moments of the proton, deuteron and neutron in 
the resonance region. The kinematic weighting of higher 
moments suppresses contributions from the low x region, 
thus minimizing the systematic uncertainty of this un- 
measured piece of the integral. Our analysis indicates 
a large TMC contamination of the CN moment /(Q^), 
which has often been used to determine the magnitude 
of higher twist effects. The alternate Nachtmann moment 
6,2 minimizes TMCs and allows a clean extraction of the 
twist-3 reduced matrix element. The previously undeter- 
mined c?2 provides unambiguous evidence of dynamical 
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Proton 
vaiue otat ^yst 


Deuteron 
value otat oyst 


Neutron 
value lotai 


Non-Singlet 
value lotai 


^ 2 
T^m 
i 2 
— m 

2 


-0.0132 • • • 0.0004 
-0.0138 0.0007 0.0009 
0.0126 0.0008 0.0020 


-0.0219 • ■ ■ 0.0009 
-0.0107 0.0012 0.0009 
0.0129 0.0016 0.0018 


-0.0105 0.0005 
0.0022 0.0020 
0.0013 0.0034 


-0.0027 0.0007 
-0.0160 0.0028 
0.0113 0.0051 


Ar2 


-0.0006 0.0008 0.0021 


-0.0090 0.0016 0.0020 


-0.0092 0.0035 


0.0086 0.0051 



TABLE III. The higher twist Ar2 = Fa - Tj evaluated at (Q^) = 1.28 GeV^ 



twist-3 effects, which arise from quark-gluon interactions. 
The first moment of §2 indicates an isospin dependence of 
HT effects in the region of small Bjorkcn scaling variable, 
assuming that the BC sum rule holds. After nearly thirty 
years of spin-dependent investigations, the 52 structure 
function remains largely unknown. But the data reported 
here and future analyses [s^ will shed light on this fun- 
damental nucleon observable. 
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